Phase dynamics after connection of two separate Bose-Einstein condensates 
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We study the dynamics of the relative phase following the connection of two independently formed 
Bose-Einstein condensates. Dissipation is assumed to be due to the creation of quasiparticles induced 
by a fluctuating condensate particle number. The coherence between different values of the phase, 
which is characteristic of the initial Fock state, is quickly lost after the net exchange of a few atoms 
has taken place. This process effectively measures the phase and marks the onset of a semiclassical 
regime in which the system undergoes Bloch oscillations around the initial particle number. These 
fast oscillations excite quasiparticles within each condensate and the system relaxes at a longer time 
scale until it displays low-energy, damped Josephson plasma oscillations, eventually coming to a 
halt when the equilibrium configuration is finally reached. 
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The way in which two Bose condensates, originally sep- 
arately formed, behave when they are connected poses 
a fundamental problem in non-equilibrium macroscopic 
quantum mechanics. Two separate condensates 1 and 
2 can be described by a relative number eigenstate \n) 



which in the phase representation may be written 
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where \(p) is an eigenstate of the relative phase 
if ^ ipi — ip2, with normalization {tfl^p') = S{ip — tp'), and 
n — {Ni — N2)/2 is the number of transferred atoms. 
The Fock state (|l|) may be viewed as resulting from a co- 
herent superposition of all possible values of the relative 
phase. This poses no problem as long as the condensates 
are physically separated, e.g. by a barrier so high as to be 
effectively impenetrable, since then no physical observ- 
able depends on if. The situation changes when the two 
condensates are brought into connection, e.g. by lower- 
ing the barrier to the point where Josephson tunnelling 
becomes appreciable. Then a nonzero supercurrent de- 
velops which depends on the phase through the relation 



n — I{ip) — {E J / h) svaip , 



(2) 



where Ej is the Josephson coupling energy. If we inter- 
pret the Josephson relation as an operator identity, 
we conclude that a system initially prepared in the state 
(|l|) will display a coherent superposition of states with 
macroscopically different values of the Josephson current 
(^. After some time, these different values of the cur- 
rent will give rise to macroscopically different values of 
the relative particle number 



n{t) — na + 



Idt. 



(3) 



The resulting macroscopic coherence is fragile and likely 
to be quickly lost by a tiny perturbation from the envi- 
ronment. After decoherence sets in, the relative phase 



is effectively measured and the system is expected to be 
describable by a reduced density matrix that is asymp- 
totically diagonal in the phase representation. We will 
show that such a density matrix is a uniform mixture of 
pure state density matrices, each of them lying in the 
semiclassical regime where phase and number may be 
taken as simultaneously well-defined. In general the sys- 
tem finds itself very far from equilibrium and, due to the 
interaction with the dissipative environment, it eventu- 
ally relaxes to the equilibrium configuration. 

The purpose of this work is to analyze such a dynami- 
cal process. We assume that dissipation is caused by the 
spontaneous creation of quasiparticles within each well 
due to a fluctuating number of particles in the conden- 
sate that is caused in turn by the coherent exchange of 
atoms. Such an exchange can, of course, occur only in 
the presence of a physical coupling, e.g. of the Josephson 
type, between the two condensates. When the Josephson 
coupling is operating, the Hamiltonian of the system at 
low energies may be written [|| ^ 



H = Ej{l - cos ip) + Ecn^/2, 



(4) 



where Ec represents (up to a constant) the effects of the 
interatomic interaction for small values of n. When Ej 
becomes nonzero, n is no longer a constant of motion, as 
explicitly indicated by the equation of motion (H) , and its 
dynamic fluctuations will induce excitation of quasipar- 
ticles within each well. The description of such a process 
will require going beyond the pendulum Hamiltonian (^) 
by including the dissipative coupling to the quasiparticle 
field. In practice, there may be other dissipative mech- 
anisms, such as incoherent exchange of single atoms Q] 
or coupling to the walls of a cavity |^ , but quasiparticle 
creation is an intrinsic mechanism which lends itself to 
an approximate analytical treatment. 

The problem which we address in this article is fun- 
damental in several respects. First, it presents an im- 
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portant case of (phase) quantum measurement that can 
be treated in detail. Second, the connection between 
independent EEC's has the essential features of a macro- 
scopic interference experiment 0, |[ ^, |l^, ^ |l^ , since 
it involves an initially random relative phase which be- 
comes well defined, in this case due to the introduction of 
a Josephson coupling with the help of the quasiparticle 
field. Finally, the establishment of phase coherence be- 
tween two recently connected condensates has bearing on 
the "Kibble problem" of formation of cosmological sin- 
gularities in the early stages of the Universe jl3j . Recent 
work on condensate coherence in optical lattices |Q sug- 
gests that the study of phase dynamics after nonadiabatic 
switching of the Josephson coupling is within experimen- 
tal reach. 

The analytical treatment which we present in this arti- 
cle requires the introduction of some approximations. We 
have already said that dissipation is assumed to be caused 
by the creation of intrawell quasiparticles caused by a 
fluctuating condensate number. Following ref. we 
take box normalization for the one-atom wave functions 
within each well. This is not an important limitation in 
the context of laser confined alkali condensates. The ini- 
tial relative particle number is assumed to be sufficiently 
large for the semiclassical approximation to be applicable 
just after decoherence takes over between distant phase 
values. We also assume that the two independently pre- 
pared condensates are identical and neglect quasiparticle 
tunneling. For conceptual simplicity, we will focus here 
on the case in which macroscopic quantum self-trapping 
takes place [|l6[ 0. This requires EcUq > AEj. In the 
future it will be desirable to extend the present analysis 
by relaxing some of these approximations. 

The picture that emerges from our theoretical study 
is schematically depicted in (Fig. |l|). There are at least 
three time scales in the connection process. First, deco- 
herence between different phase values takes place after a 
few atoms have been exchanged. This process effectively 
measures the macroscopic phase. The second time scale 
is the Bloch oscillation period around the initial value 
of the relative particle number. These oscillations are a 
manifestation of the ac Josephson effect due to a chemi- 
cal potential difference. Finally, the system relaxes in a 
longer time scale towards the equilibrium configuration. 
Damping is due to the irreversible energy transfer from 
the quickly oscillating number-phase degree of freedom 
to the quasiparticle field within each well. 

First we perform a semiclassical study in which the 
relative phase and number are taken as c-numbers while 
the quasiparticle field is treated quantum mechanically. 
The semiclassical assumption is later justified from a full 
quantum treatment which explicitly shows how the quasi- 
particle field transforms the reduced density matrix of the 
phase-number coordinate from that of a pure Fock state 
to a mixture of Gaussian wave packets. This is the pro- 
cess that effectively "measures" the phase and makes it 



well-defined in a particular experimental realization. 

At the semiclassical level (where the question of iden- 
tity or not between the states |(/?) and \lp + 2ti) can be ne- 
glected), it is clear from Eq. (^), that the problem maps 
on to that of the the "Bloch oscillations" of a particle 
in a tight-binding lattice which is subject to a harmonic 
potential, with n playing the role of (discrete) coordi- 
nate and that a of a wave vector. A simple semiclas- 
sical analysis reveals that a particle placed with initial 
"wave vector" Lp^ at a high point no of the harmonic 
tight-binding lattice, undergoes periodic motion. Specif- 
ically, for EcUq Ej, and with an appropriate choice of 
time origin, we may write 



n{t) = n + p sui{u!ot) , 



(5) 



where p = (Ej/nEc), luq = Ecfi/h, and, initially, n ~ ng. 
Equation ^ describes Bloch oscillations undergone by a 
quantum particle in a tilted lattice. The resulting self- 
trapping is not stable because the rapidly oscillating n{t) 
induces excitations in the quasiparticle field, which acts 
as a dissipate environment for the macroscopic degree of 
freedom |l^ . These oscillations are very fast and one 
is rather interested in the slower dynamics of the time 
average fi as it experiences damping. 

Since quasiparticles within each well are sensitive to 
the number of particles in the condensate to the extent 
that they are but linear deviations from the mean field de- 
scription 1 20 , it is clear that (^) yields a time-dependent 
perturbation that will excite quasiparticles. The time- 
dependent Hamiltonian in e.g. the left well is |^ 
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where g is the coupling constant, N is the total particle 
number, and V is the volume of one box. The Hamilto- 
nian for the right well looks like ^ with n(t) replaced 
by by —n{t). Eq. (0) may be rewritten as 



H{t) = Ho + Apsm{ujt), 



(8) 



where Hq is the unperturbed "bath" Hamiltonian [with 
n{t) replaced by n] and 
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A = ^^^ ]{2clck + 4.cL^ + c^kCk)- 



(9) 



By applying linear response theory with n treated as 
constant, we obtain for the dissipated power within one 
well: 



where 
X"(^) 



W - (pV2)^ox"('^o), 



^ I dte^-'{[Ait),AiO)]) 



(10) 
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In the third line, we have taken the low frequency ap- 
proximation huj <^ mc^, c being the sound velocity and 
771 the atom mass pl[ . Here Uk and Vk are the coherence 
factors of the Bogoliubov transformation, so that, in this 
limit Uk - Vk ^ {nVck/Ngf/'^, where N/2V is the 
atom density in one well. We note that x"(w) goes like 
(jj'* at low temperatures and like w'^ at high temperatures. 
The second line of Eq. ( pi] ) indicates that the quasipar- 
ticles are excited in pairs, which can be ascribed to the 
pairing structure of the depletion cloud |3[ ^ . 

The energy that is dissipated into the quasiparticle 
field has to be taken from the energy of the macroscopic 
pendulum (^ which, under the self-trapping assumption, 
is dominated by the interaction term. Thus, we write 



2W = -Erfin, 



(12) 



where the factor of 2 accounts for the presence of two 
wells. We derive the differential equation 



= -Cn coth 



Wbt) ' 



(13) 



where C = 2TTr]E'j/nN'^Ec, with Ec = 2g/V = 
STrTi^a/mV and rj = {No? /2'kVY^'^ the small parame- 
ter of the dilute gas. At low temperatures {ksT <C Ecfi) 
we get 



nit) 



no 



1 + 7ri' 

while, at high temperatures (ksT 3> Ecfi), 
n{t) — no exp(— 7-rt). 



(14) 



(15) 



In (|4|) and (|T5|), 7r is, respectively, the low and high 
temperature limit of the general expression 



hC ( hujo 

7. = -^0.0 coth 



(16) 



where wq refers to no- 

As the dissipation process continues, the energy of the 
macroscopic pendulum decreases until it eventually un- 
dergoes low energy harmonic oscillations. In this limit, 
the analysis of the dissipation process is similar to that 
described above, with t^o replaced by the Josephson 
plasma frequency wjp — \/EjE^/h and with the pendu- 
lum energy written as n^Ec/2, where n is the amplitude 
of the number oscillations. Since lu jp does not depend 
on 77. for n N, we obtain for both low and high tem- 
peratures 



h{t) = 77(0) exp(— 7rt), 



(17) 



7r being given by (16) with wq replaced by ujjp. 

It remains to justify the semiclassical approximation. 
For that we have performed a full quantum calculation of 
the phase dynamics right after the connection has been 
established. The direct reversal process 77, —n may 
be safely neglected. We generalize quantum dissipation 
techniques |2^ to the case of quadratic coupling in the 
bath coordinates [see Eq. (||)]. After some lengthy al- 
gebra, we find that, in the phase representation, the re- 
duced density matrix simply evolves as 



sin-^ 



ip-tp 



where 



7d = 7r"o coth 



2kBT'' 



(18) 



(19) 



For 7rft » 1, Eq. ( [l8| ) becomes effectively equivalent to 

\p{v,v';t)\ - {2tt)-^ eM-{v - V' flat], (20) 

and hence 7^ is identified as the decoherence rate. It is 
possible to prove that the reduced density matrix ( [20| ) is 
a mixture of pure state Gaussians ~ exp[— 2((/? — 0)'^^cLt] 
with 9 uniformly distributed between and 27r. Thus we 
assert that the density matrix ( pO| ) describes a quantum 
system in which the variable (p is measured with preci- 
sion ~ [idt)^^'"^ with the average result chosen among a 
menu of uniformly distributed possible values. Here the 
role of the "measurement apparatus" is played by the 
quasiparticle field with its many degrees of freedom. 

At low temperatures, we obtain ^d/lr = i^q. Now 
we note that 27^ is the relaxation of the total energy 
Ecn^/2 and that E^no is the variation in energy due to 
the exchange of a single atom. We conclude that the 
ratio between the decoherence and the energy relaxation 
rates is the ratio between the total energy and the energy 
variation associated to the transfer of one atom. Thus 
decoherence is achieved as soon as \n{t)\ decreases by a 
few units (H] . Once this has happened, Eq. (|20| ) apphes 
and the system may be treated semiclassically with an 
initial value of the phase that is chosen at random and 
which is defined with increasing precision as time evolves 
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In other words, the decoherence rate can be 
identified with the rate of net atom exchange. On the 
contrary, an appreciable energy relaxation requires the 
exchange of many atoms. This picture is consistent with 
the idea that the net exchange of a single atom plays the 
role of the excitation of one quantum of the dissipative 
bath 1^. 

Once we can describe our system as a mixture of semi- 
classical wave packets, linear evolution applies to these 
even in the presence of dissipation. Since the transla- 
tional motion of a Gaussian wave packet is semiclassical, 
we apply the method initially described with the confi- 
dence that it has been justified by a full quantum treat- 
ment. 
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The final picture of the connection process is one which 
involves several time scales. Decoherence in the phase 
representation sets in after the exchange of a few atoms 
has been completed. For a given experimental run this 
can be interpreted as the measurement of the phase in 
which a given value is chosen with an accuracy that im- 
proves with time. After the phase is defined, a semiclas- 
sical treatment suffices to conclude that the macroscopic 
degree of freedom will undergo very fast Bloch oscilla- 
tions that will dissipate energy into the many degrees of 
freedom of the quasiparticle field. Ultimately this Bloch 
oscillations become low energy Josephson plasma oscilla- 
tions that are also damped by the quasiparticle field. At 
long times we expect that the system density matrix will 
be the thermal equilibrium one which corresponds to a 
total energy of Ecn'^/'2,. For large V the corresponding 
temperature is much smaller than the critical tempera- 
ture. 

The main message is that the depletion cloud, de- 
scribed by the quasiparticle field, can be an intrinsic 
agent for the decoherence and measurement of the macro- 
scopic phase. We have studied the most fundamental 
intrinsic mechanism. Incoherent exchange of thermal 
atoms Q and coupling to the walls of a cavity will 
yield additional contributions to both the decoherence 
and the energy relaxation rates. 
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FIG. 1: Schematics of the connection process in the number 
representation: After a time t ~ l/ld a semiclassical wave 
packet forms that undergoes Bloch oscillations thus exciting 
quasiparticles. For t > l/7r, the system relaxes its energy, 
eventually performing underdamped Josephson plasma oscil- 
lations. 



